Entanglement of bosons in optical lattices 
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Entanglement is a fundamental resource for quantum information processing 1 1 2 1 which occurs naturally in 
many-body systems at low temperatures. The presence of entanglement and, in particular, its scaling with the 
size of system partitions |3 - 5| underlies the complexity of quantum many -body states |5-7|. The quantitative 
estimation of entanglement in many-body systems represents a major challenge as it is held to require either 
full state tomography |8 |, which scales exponentially in the system size, or the assumption of unverified system 
characteristics such as its Hamiltonian or its temperature. We adopt recently developed approaches for the 
determination of rigorous lower entanglement bounds from readily accessible measurements |9 10 1 and apply 
them in an experiment of ultracold interacting bosons in optical lattices of approximately 10^ lattice sites. We 
use this approach to study the behaviour of spatial entanglement between the sites when crossing the superfluid 
to Mott insulator transition and when varying the temperature. This constitutes the first rigorous experimental 
large-scale entanglement quantification in a scalable quantum simulator. 



Background - The last few years have seen experiments to- 
wards the verification of the presence of entanglement in a 
variety of physical realizations of many-body systems. Multi- 
particle spin entanglement of distinguishable particles has 
been created and studied experimentally for up to 14 sites 
in ion traps |11-13| and up to eight sites in photonic setups 
|[T4l by means of entanglement witnesses that determine the 
presence of entanglement. In ultracold neutral atomic gases, 
entanglement between indistinguishable particles with two in- 
ternal degrees of freedom was generated by squeezing of the 
total (pseudo) spin and its presence verified by spin- squeezing 
inequalities |15-18|. For ultra-cold bosons with two internal 
degrees of freedom in optical lattices, entanglement between 
lattice sites was created by controlled collisions and qualita- 
tive evidence for its presence was found 1 19|. Experimental 
evidence for the entanglement between the spins in magnetic 
materials was given by comparison of the neutron scattering 
structure factor to a classical description |20|. 

A key challenge that remains to be addressed however 
concerns the quantitative determination of upper and lower 
bounds on the entanglement that is present in the system with- 
out resorting to assumptions concerning unverified system 
characteristics such as its Hamiltonian or its temperature. This 
can be achieved by adopting a simple but powerful principle 
f9", '2r-'23 1 in combination with methods from optimization 
theory: Given a set of observables, we consider all density 
matrices that are compatible with experimentally measured 
expectation values of these observables. Amongst these den- 
sity matrices, we find the one with the least (largest) amount 
of entanglement as quantified by a suitable entanglement mea- 
sure 1 1 1. In this way we determine a lower (upper) bound on 
the entanglement that must have been present in the state that 
gave rise to the observed expectation values. These bounds do 
not require any other assumptions on the system. 

Fundamentals - The multi-partite entanglement that we ob- 
serve and quantify in our experiment arises in a periodic op- 
tical potential (optical lattice) that hosts massive bosonic par- 




FIG. 1 : We investigate entanglement in an optical lattice filled by 
ultracold bosons. At large lattice height, the ground state is a pure 
product state and no quantum correlations exist between any subset 
of sites (A, B, C, ...). Such a state can be created by only performing 
physical local operations (those that respect the mass superselection 
rule locally, supLO operations) and allowing the parties associated 
with the subsets to communicate classically (CC). Decreasing the lat- 
tice depth, the resulting state may not be created by such operations 
anymore - the state becomes a resource of value 8 with which the 
parties may overcome their locality restrictions, i.e., entanglement is 
created. 



tides (see Fig. [T]). Hence, the non-local correlations are un- 
avoidably intertwined with the superselection rules that pro- 
hibit the formation of coherent superpositions with different 
particle numbers. As a consequence, before we can provide 
a rigorous theoretical and experimental quantification of the 
entanglement in the system, we need to clarify its nature. We 
achieve this from the viewpoint of entanglement as a resource 
LU. 

Physical constraints, fundamental or practical, impose lim- 
itations on accessible physical operations (MlES |. For exam- 
ple, the locality constraint that expresses the inability to ex- 
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FIG. 2: The minimal entanglement E{k) consistent with the measured momentum distributions of an atomic gas in a three-dimensional optical 
lattice tuned across the quantum phase transition from superfluid to Mott insulator (lattice depth s = 9, 12, 15, 18, 21). For each value of the 
momentum k in the first Brillouin zone, E{k) provides a lower bound to the entanglement present in the bosonic sample. The total number of 
particles was approximately 3.5 x 10^ for each measurement. See Fig. 2 for total particle numbers. 



change quantum particles between distant laboratories A and 
B prevents the execution of quantum gates between A and B. 
Such constraints in turn imply the existence of resources — 
here entangled states — that, when consumed, allow for the 
realization of operations that are impossible under the given 
physical constraints [26 ] . The use of entanglement in telepor- 
tation for example allows for the realization of general quan- 
tum gates between A and B ||27]|28l. This aspect captures the 
multi-partite nature and non-local correlations in the system. 

For massive bosonic particles an even more fundamental 
constraint concerns superselection rules for massive indistin- 
guishable particles: Physical operations cannot create coher- 
ent superpositions of different particle numbers. In connection 
with the locality constraint this requires that all physically al- 
lowed local operations must commute with the local particle 
number operator. 

The fact that both constraints, locality and super- selection 
rules, need to be considered simultaneously for massive 
bosonic particles leads to a refined picture of non-local re- 
sources and thus entanglement |29-32|. In a multi-partite 
system, as in Fig. [T] in which two or more parties aim to ex- 
change quantum information, but are restricted (i) to only act 
locally on their respective quantum system and communicate 
classically (LOCC), and (ii) to perform only operations pre- 
serving the local particle number operator, they may (in the 
two-partite case) only prepare states g of the form 

B = Y.P-e^A^Q^B\ (1) 

n 

(n) (n) 

where g\ ' (^^ are density operators of subsystem A (B), 
{Pn} is a probability distribution and the local states must 
commute with the local particle number operators Na and 
Nb. respectively, i.e., = = 0. All 

such states are here collected in the set S. All other states, 
so states not in 5, become a resource, to be used to overcome 
the constraint imposed by locality and/or superselection rules 
(see Appendix [A| for an illustrative example). 

Entanglement Quantification: Theory - Having clarified 
these fundamental issues, the key point is now that we turn 
this qualitative appreciation of the value of these states into 
a mathematically and physically well defined quantifier that 



may then be determined experimentally. It is a crucial require- 
ment that the value of this quantifier does not increase on aver- 
age under LOCC operations that respect local superselection 
rules (supLOCC in short). Such a quantifier is then denoted a 
supLOCC monotone. To this end, we define W as the set of 
entanglement witnesses W [33] satisfying (i) tr[Vr^] > for 
all ^ G 5 and (ii) the operator inequality W ^ N > 0. Then 

f(^) = max|o,- inf tr[Wg]} (2) 

is a supLOCC monotone (see Appendix[B]for a proof and note 
the similarity to entanglement monotones expressed as opti- 
mization over witnesses for spin systems in Ref. |34|). Note 
that, for any state g, its entanglement S{g) is upper bounded 
by the mean total number of particles (TV) = tv[Ng], pro- 
viding a figure of merit for the lower bounds we will present 
below. This definition may be straightforwardly extended to 
more than two parties and we now set out to quantify the en- 
tanglement under superselection rules contained in states of 
bosons in optical lattices when each lattice site constitutes a 
party, i.e., labelling the sites of the three-dimensional lattice 
by i = (n, ^2, ^3), states in S are of the form pn g^l^^ 
with [g^l^\ni] = 0, where hi is the number operator for the 
lattice site i. In this way, we quantify the entanglement shared 
between sites of the lattice (as opposed to between particles). 

Needless to say, Eq. ^ is exceedingly hard to compute an- 
alytically or even numerically, especially in the many-particle 
system (where each constituent is in addition being described 
in an infinite-dimensional Hilbert space) that we consider 
here. However, rather than aiming for exact values we will 
follow Refs. (21 ED [22 [3l to derive lower bounds to Eq. ^ 
which, after the introduction of essential aspects of the experi- 
mental setup, we will demonstrate to be obtained from readily 
accessible measurements. 

Experiment - In the experiment we quantify the entangle- 
ment of a system of ultracold interacting bosons in a three 
dimensional lattice potential. An almost pure 13511 Bose- 
Einstein condensate of ^ 3.5 x 10^ atoms of ^^Rb is prepared 
evaporating a sample of atoms in the \F = l^Mp = — 1) 
state in a hybrid trap composed of a focused red detuned laser 
beam (optical dipole trap ODT) propagating in the horizontal 
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plane plus a quadmpole magnetic field. As first demonstrated 
in 1 36] , when the focused laser beam is slightly offset ver- 
tically from the center of the quadrupole magnetic field, the 
atoms experience a harmonic potential with cylindrical sym- 
metry. In our system the resulting frequencies are ^50 Hz 
and ^8 Hz in the radial and axial directions respectively. An 
optical lattice (OL) potential with lattice constant a = A/ 2, 
generated by three counter-propagating red-detuned beams 
(with wavelength A = 830.3 nm and waists ^ 180 fim) 
is slowly superimposed to the sample by performing an ex- 
ponential ramp in tramp = 140 ms. The final amplitude of 
the lattices Vql = s£^r (where = /(2m\^) is the re- 
coil energy associated to the absorption of a lattice photon by 
an atom with mass m) has been calibrated with an accuracy 
of ±10% through lattice amplitude modulation spectroscopy 
1 37 1 and can be varied from 5 = to s = 30. In this way, 
we realize a many-body state of bosons in a three dimen- 
sional cubic lattice (subject to the harmonic trapping poten- 
tial), the entanglement of which we are interested in quan- 
tifying. The system's Hamiltonian is well approximated by 
the Bose-Hubbard Hamiltonian [38 1 (note that this informa- 
tion/approximation will not enter into our quantification of the 
entanglement). At sufficiently low temperatures, when the lat- 
tice depth is s > 15, the ratio between the interaction energy 
U of two atoms in the same lattice site and the tunnel energy 
J between two adjacent lattice sites is large enough to obtain 
a Mott insulator as firstly demonstrated in 1391. Due to the 
overall harmonic confinement an inhomogeneous Mott insu- 
lator is obtained with regions with different filling. At s = 27, 
we estimate ~70% of the atoms (in the outer region of the 
sample) to be in a Mott shell with 1 atom per site and ~30% 
(in the central region) to be in a Mott shell with 2 atoms per 
site. 

We now describe the experimental procedure — time-of- 
flight measurements — which will give access to an observable 
allowing us to lower bound the entanglement of the bosons in 
the optical lattice. After a holding time thoid = 5 ms, we si- 
multaneously switch off both the trap and the optical lattice. 
The cloud then expands freely for a time ttof = 21 ms be- 
fore we measure the column density through absorption imag- 
ing. In practice, we measure the optical density fi{x,y) of the 
atomic sample integrated along the direction of the imaging 
laser beam, and we extract the real atomic (column) density 

n{x, y) =a (/i(x, y) - /io) (3) 

with a being a pre-factor related to the imaging calibration 
and the effective size of the CCD square pixel, as discussed 
in Appendix |D] and /io is the background noise level of the 
image, that is mainly due to residual fluctuations in the laser 
intensity of the imaging beam. 

Entanglement Quantification: Experiment - From the col- 
umn density n{x,y) we can extract a lower bound to the 
entanglement £ as follows. The column density at position 
{x^y) = r = htk/m (where k = {k^ ky) is the quasi- 
momentum in the lattice) after a time of flight t is given by 
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FIG. 3: Minimal entanglement E(k) consistent with time-of-flight 
measurements as a function of the lattice depth s. The reported data 
correspond to the average of the lower bound E(k) over 5x5 pixels 
centered at k=(7v/a, n/a). Lines are guides to the eye. Error bars 
include statistical and systematic errors due to the uncertainty in the 
lattice depth and in the image acquisition and processing, see Ap- 
pendix |D] for details. Circles in grey show the corresponding mean 
total number of bosons (without error bars for clarity, relative error 
^ 8%) which upper bounds the entanglement as discussed in the text 
after Eq. 

iolEIl 

h{k) = f(k) blbje^'^^^-^^e'^^^^'-^'l (4) 

where P- (hi) is the creation (annihilation) operator of a par- 
ticle at the lattice site i and f{k) > may be obtained from 
a numerical band-structure calculation (See Appendix [C] for 
details). 

It is straightforward to show ^ that W{k) := 
n{k)/f{k) — TV G W, where TV is the total particle number 
operator, so that it is a witness as required in Eq. ([2]). Then, 
for any state g, we have a lower bound to its entanglement 
content in terms of (TV) = tr[^TV] and {n{k)) = tv[gn{k)]\ 

£{e) > max{0, (iV) - =: E{k), (5) 

which holds for all k. Note that there are no assumptions: The 
entanglement of any state is bounded from below by Eq. ([5]). 
The Hamiltonian governing the system, the temperature, de- 
tails of external potentials, or even the system being in equi- 
librium, need not to be assumed. 

As E(k) is a lower bound to the entanglement £{q) for 
all k, averages over an area A, jj^dkE{k)/\A\, also provide 
lower bounds. We use this fact to account for the finite reso- 
lution of the camera and to incorporate symmetries (see Ap- 
pendix |D] details). For ease of notation, we denote this lower 
bound also hy E{k). 

We determine the minimal entanglement E{k) for different 
optical lattice depths across the superfluid to Mott insulator 
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FIG. 4: Lower bound E(Jk) as a function of temperature for three 
different values of the lattice depth [s = 6 (filled boxes), s — 12 
(filled circles), s = 18 (filled triangles)] following the caption of 
Fig.|3] Corresponding non-filled symbols in grey show the mean total 
number of bosons. Horizontal error bars indicate the uncertainty due 
to the calibration of the imaging system (error bars on mean total 
number of atoms not shown for clarity, relative error was ^ 8%). 



transition. Before presenting our experimental results, let us 
consider the two extreme cases analytically. For ultra-deep 
lattices {s oo) and at zero temperature, the system will 
be in a Fock state |nin2 . . . ), for which E = E(k) = by 
definition. For very shallow lattices, when tunnelling domi- 
nates over the on-site repulsion and one neglects the latter, the 
ground state with N particles of the translationally invariant 
Bose-Hubbard Hamiltonian is proportional to 6t)^|vac) 
and one finds E{k) ^ N at k = k := (7r/a,7r/a). We thus 
expect the entanglement to decrease when increasing the lat- 
tice depth. 

Results - In Fig. [2] we show E{k) in the first Brillouin zone 
for s = 9, 12, 15, 18, 21. For each value of s we collected 
^ 40 absorption images in order to reduce the statistical error 
on the determination of the entanglement. Relative shot-to- 
shot spread of the atom number is lower than 10%. The lower 
bound of the entanglement decreases as the system crosses 
the transition from the superfluid (lower s values) to the Mott 
insulator phase (s > 15). This behavior can be seen better 
in Fig. |3] reporting the value of E{kx^ky) averaged over a 
box of 5 X 5 pixels around k, where we expect [9J and found 
the bound to be largest. Details of our imaging system are 
explained in Appendix [D] where we also describe the error 
analysis. 

As the entanglement of the system is expected to decrease 
with increasing temperature f9|, we also perform measure- 
ments fixing the optical lattice depth s and varying the tem- 
perature of the atomic sample. The determination of the tem- 
perature inside the lattice is still challenging [42 , 43] while its 
measurement in a harmonic potential, i.e., before raising the 
lattice, is routinely done. Here, we refer to the temperature T 
before the loading of the optical lattice. In practice, in order to 



realize samples of different temperatures, we perform optical 
evaporation in the ODT to different values of the power Pq dt 
and then we increase the power of the optical dipole potential 
up to a fixed value Pi. This procedure allowed us to obtain 
temperatures from 40 nK up to 100 nK in the same final har- 
monic potential before loading the lattices (see Appendix [E| 
for details). In Fig.|4]we show the behaviour of the minimal 
entanglement for different temperatures and for three differ- 
ent values of the lattice depth corresponding to a superfluid 
ground state (s = 6), a Mott insulator phase (s = 18), and the 
crossover region (s = 12). As expected |[9|, with increasing 
the temperature, the minimal entanglement consistent with the 
measurements decreases. 

Conclusions & Outlook - We have quantified experimen- 
tally the multi-partite entanglement of a system of interact- 
ing bosons in an optical lattice through routinely done mea- 
surement of the atomic density profile after expansion. As 
the Hamiltonian and the ensuing dynamics of such a system 
can be controlled, it constitutes a bosonic quantum simulator, 
naturally supplying the resource entanglement at low temper- 
atures. Our estimation of the entanglement is rigorous and 
without unspoken assumptions and provides a quantitative in- 
sight into the structure of the many-body state. In essence, 
we have answered the question, "Which is the least amount 
of entanglement that is consistent with given measurements?" 
["21 1 . The strategy we implemented for this estimation is suf- 
ficiently general to allow for its adaption and application in 
a wide variety of experimental settings that arise naturally in 
quantum science. Indeed, this principle may also be general- 
ized to other quantities. One may for example ask what is the 
maximal entropy consistent with given measurement results 
and by answering this question place rigorous, assumption- 
free, upper bounds on the entropy of a quantum many-body 
systems. Of course, our approach is not restricted to bosonic 
systems but may also be applied to fermionic and spin systems 
or mixtures of bosonic and fermionic atoms, thus providing 
quantitative information about complex states of matter. 
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Appendix A: An illustrative example 

We illustrate the resource character of states that are not of 
the form Eq. ([T]) by means of an example taken from Ref. 13Q|. 
Suppose a single classical bit is encoded in the relative phase 
of the two states 

I±) = ^(|0>a|1)b±|1)a|0)b). (Al) 

Remarkably, if Alice and Bob are constrained by LOCC and 
local particle number conservation, they are unable to distin- 
guish these two states - the bit ± is hidden from them. They 
may learn the bit, however, in one of two different ways: (i) if 
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they share entanglement in the form 

1 ^ 

which has (V^|iV|^) = N, trB[|V^)(V^|] and trA[|V^)(V^|] com- 
mute with the respective local number operators, but it may 
not be written as in Eq. ([T]). Or (ii) if they share the state 

(z = \a\&^) 

g(x dcl)\z)A{z\0\z)B{z\, |2;)(xV^|n), (A3) 
Jo ^0 

which is of the form in Eq. ([T]), commutes with the total num- 
ber operator, and has (TV) = tr[7V^] = 2|ap. 

Now, the success probability of learning the phase of |±) 
increases with (TV) (approaching unity as (TV) oo), or, in 
other words, the value of these resource states increases as 
their mean total number of particles with (TV) increases. For 
Itp) the entanglement increases while for g the local violation 
of particle number conservation increases. 

Let us evaluate the lower bound in Eq. ([5]) for these states. 
Suppose that the lattice consists of the two sites a = and 
6 = (1 0)\ which we associate with Alice and Bob, respec- 
tively. For the states in Eq. ( |A2| ) and Eq. ( |A3| ), we find 



E 



2\a\'=tr[Ngi 



(A4) 



respectively. Hence, for both states, we have a lower bound to 
the entanglement that is increasing in (TV) (thus capturing the 
value of these states for the above data-hiding protocol) and 
for the second state the bound is in fact exact as f (^) < tr[TV^] 
for all states g. 

Appendix B: A monotone under SSR-LOCC operations 



Here we show that 



f(^)=max|o,- inf tr[V^^]| 



(Bl) 



is a monotone under LOCC operations commuting with the 
local particle number operators. Let the Hilbert space H be 
a direct product of Q parties, H = ®^=i^s- LOCC op- 
erations with respect to this partition are operations taking 
density matrices g to J2kPkQk^ where g^ = AkgA\/pk, 
Pk = tr[A/e^A^] and the are of the form = (8)s=i 
and fulfil AlAk < 1 and [A^, fis] = 0. £{g) is an entan- 
glement monotone if 



^PkE{gk)<£{g). 



(B2) 



For all ly G W, we have tr[Ty^/e] > -tr[Ngk] > -oo. 
Hence, the infimum exists and we denote it by . Now let 



e > 0. Then Ej. -\- e is not an infimum and therefore there is 
an Wk,e ^ W such that tr[Wk,eQk] < Ej, ^ e, i.e.. 



^^PkEigk) = - PkEk < ^ Pfe(e-tr[iyfe,e^fe]) 

k 

Ek<Q 

( E ^lWkAk)Q . 



k 

Ek<0 



< e - tr 



k: Ek<0 



which is upper bounded by 



y^PkEigk) <e- inf tr[W'^], 

k wew 



(B3) 



(B4) 



as the operator in brackets in Eq. (B3 ) is a member of W: For 
all k and all e, we have Wk^e + ^ > 0, i.e.. 



0< 4(^fe,e+^')ife< ^{WkAk^N 



k 

Ek<0 



k 

Ek<0 



(B5) 



as [i/e, TV] = for all k. Now let a eS. Then 

{Y ^iWk,eAk)a = Y tr[Wk,eAkaAl], (B6) 



tr 



k 

Ek<0 



k 

Ek<0 
it 



where, up to normalization, AkaAl G 5, i.e., the above is 
lower bounded by zero (as Wk^e ^ W) and we hence have 
that for all e > 



YPk^(^k) < e - mf ir[Wg] < e + E{g), 



(B7) 



which implies that £ is an entanglement monotone. 



Appendix C: Density after time-of-flight 

We set out to derive an expression of the atomic column 
density after free evolution, i.e., after evolution under the 
Hamiltonian 



H 



2m 



(CI) 



In order to connect the atomic column density to observables 
in the lattice, i.e., before the free expansion, we expand the 
field operators in Wannier functions of the lattice 



(C2) 



Here i is a multi-index containing the lattice site and the band 
index. The density operator after evolution under H for a time 
t reads 

= ^<(r)^,•(r)e^^^/^6l6,•e-^^^/^ ^^^^ 
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Now, due to the lattice geometry, the Wannier functions fac- where we denoted {f^g) = / dr f''{r)g{r). The Hamilto- 
torize and are the same for each spatial direction, Wi{r) = nian is diagonal in the basis of the a^, 
Wi^{x)wiy{y)wi^{z). Owing to orthonormality, we hence 
find for the column-density operator n(x, t) := J dzn(r,t) 
after time-of-flight t 



n{x,y,t)= wl{x)wj^{x)w;^{y)wj^{y)bl{t)bj{^^ 



which implies 



(C6) 



To compute the time-evolution bi{t) = q^^h /h'j^^^-itH /h 
the bosonic annihilation operators, we expand in orthonormal 
and complete plane waves uj{r) = q^^^p-^^/l j ]j>/2 ^ 



dp, 



^^^^ i.e., using Eqs. ( |C5] ), 



(C7) 



Due to orthonormality, the bosonic annihilation operators are 
related as 

hi = Y{wi,ujp)ap, CLp = ^(cjp, Wi)hi, ^qi^-^ 



p J 



Hence, 



n{x,y,t)= ^ wl{x)wj^{x)w;^{y)wj^{y){up,Wi){wj,Uq)^^^^ ^^^^ 
where completeness and orthonormality imply ^ (cjp^ , Wi^ ) (k;^^ , ujq^ ) = and k;^ (x) (i^j , ujq) = ujq{x), i.e.. 



n(x,y,t)= ^ <(xH„(xV;^(yH^(y)e"r&i(p'-''')ataq. (C9) 



Now, again using Eqs. ( |C5| ), we arrive at an expression for 
the column density in terms of the hi 



i{x,y,t)= "Y gl{r,t)gj{r,t)h\hj, 



iz=jz 



(CIO) 



where gi{r, t) = gi^ (x, t)gi^ {y, t), with 

gi{x,t) = Yujq{x){ujq,Wi)e~'^^'^ 
q 



(Cll) 



and we used completeness and orthonormality to arrive at 
J2p, i^iz ^ ^Pz ) i^Pz 5 '^jz ) = ^iz dz ' Finally, considering the 
properties of the Wannier functions, and after some algebra, 
we let L ^ oo to obtain 



gi{x,t) 



2 

J — ( 



27r ^-oo 



where r = and 

1 

^n(^) = ^ / drwo,nMe-^^'^' (C12) 

V ZTT J-oo 



is the Fourier transform of the Wannier function of the n'th 



27r 



band centred at zero. Eq. ( |C1Q| ) with gi as above is the exact 
expression for the column density at {xy). All the involved 
functions may be obtained by a numerical calculation of the 
Wannier functions. 

In the experiment, only the lowest band is occupied and we 
omit the band index from now on. In the stationary phase 
approximation, which is valid for 1 <C r (~ 1.8 x 10^ in our 
experiment), one has 

f,{x)^{l-i)^w{^x-^i). (C13) 
Finally, approximating w{^x — ^i) w{^x) |41|, yields 
n{r = ^k,t) = f{k) J2 e'l'^*'(^-^')+"'(^"-^')/^l6t6,-, 

iz=jz 
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for the column density at r = {xy) = ht{kx ky)/m = 
htk/m = TO?k /{2ti'^) after a time-of- flight t. Here, 



function / by /(x, y) ^ f{xi,yj) for each (x, 7/) G (A^, A^), 
we introduce the simpHfied quantities Ea and Eij 



1 

w{^k) = I dr wo{ar)e 



-ikar 



(C14) 



and Wo is the Wannier function of the lowest band centred at 
zero. 



Appendix D: Image and error analysis 



We intend now to analyze the measurement of 



E{x,y) = {N) 



{n{x,y)) 
f{x,y) 



(Dl) 



where we now work in real-space coordinates, i.e., 
{n{x^y)) = J dz {n{x ^ y ^ z)) , where {n{x^y^z)) is the ex- 
pectation value of the density distribution of the atom cloud at 
r = {xy z) after time-of-flight and 



f{x,y) 



(D2) 



as in Eq. ( p4| ). 

Due to the spatial discretization of the CCD sensor used in 
the experiment, we define the discrete function 



E' 



1 

A2 

(TV) 



dx / dyE{x^y) 

Idxldy^^^^ 



(D3) 



where (i, j) denotes the index of each pixel, centered on 
{xi,yj),dind /^i^j = [xij - A/2,Xij + A/2]. A = 2.78 /im 
is the effective pixel size which takes into account the physical 
pixel size and the magnification of the imaging system. The 
total number of atom is given by = j riij. Note that 
the quantity E^j is still a lower bound for all (z, j). 

For each acquired image, we incorporate the symmetry of 
the observable n{k)/ f(k) by averaging over pixels corre- 
sponding to {kx^ ky), {kx ± 27r/a, /c^ ± 27r/a), and the sym- 
metry of the experimental setup by averaging also over the 
four points (i/c^^, ^ky). For Figs, [s] and |4j we additionally 
consider the average of E^j on a subset of 5 x 5 pixels (cor- 
responding to twice the width of the point spread function of 
the imaging system) centered around k G [— 27r/a, 27r/a] 
where a is the lattice spacing. If we define a set of pixels A 
over which we perform the average, the quantity 



(D4) 



is also a lower bound to the entanglement. 

Actually, we do not have access to the quantity 
n{x^y) / f{x^y) to be integrated in Eq. (D3). Approximat- 



ing (and taking the error into account below) the continuous 



E'a ~ Ea 



(D5) 



where gi^j = l/f{xi,yj) and {fiij) is the expected number 
of atoms recorded by pixel (i, j). 

In the experiment, we use a running average of about 40 
density profiles for each value of lattice depth s and tempera- 
ture T. The empirical average of Ea over a set of M images 
is also the best estimation of the entanglement bound 



M 



Ea 



X/fY.^A^ - \A\ E ^^'i' 



M 



(D6) 



n=l 



n = 1, . . . , M being the image index. We proceed by analyz- 
ing the sources of uncertainty when estimating E'j^, i.e., the 
systematic uncertainty related to ni^j and gi^j, and the statis- 
tical contribution associated to shot- to-shot variations of rii^j. 
The statistical uncertainty can be estimated as 

M 



( stat\2 
i'^E^ ) = 



M(M 



(D7) 



The approximation g{x,y) gij used in Eq. (D5) for x G 
Ai^y e Aj introduces a systematic error. We find 



< 



dx dy 

J A, JAj 



{n{x,y)) 
A2 



g{x,y) -gij 



. (D8) 



From the mean value theorem we have that 

9ij -9{x,y) = {dxg){a,b){xi - x) + {dyg){a,b){yi -y), 

where (a, 6) = (1 — c)(x, y) + c(xi, yi) for some c between 
and 1. Hence, for x G A^, G A^, we find 

-9{x,y)\ < ^ max \Vg{a,h)\ =: ^6^,^, 
aeAi,beAj ^2 

and with this result (|D8|) becomes 



< 



V2A 



(D9) 



Thus, assuming a flat error distribution, the resulting standard 
error is given by 



VGA 



(DIO) 



Now we analyze in more detail the other systematic con- 
tributions. In the following, we write i = {ij). The gi are 
obtained from the Wannier function of the optical lattice and 
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thus the error depends on the uncertainty as = 0.1s we have 
in the estimation of the lattice depth s, i.e., 

< = (DID 

(n) 

Now we discuss the systematic error on . In the ex- 
periment, we measure the optical density distribution by ab- 
sorption imaging. More specifically, we record on a CCD 
camera the intensity profile of a resonant probe laser beam 
interacting with the sample. The absorbed light intensity la, 
integrated along the z direction (the direction of propagation 
of the probe beam), as given by the Beer-Lambert law, is 

la = Iq (^1 — exp[— cr J h{r)dz]^ , with Iq being the incident 
intensity and a the resonant absorption cross section given by 
a = 3AQ/(27r), where Aq is the wavelength of the resonant 
transition. For circularly-polarized light on the transition we 
used for ^^Rb, we have a = 2.907 10"^^ m^ |i44J. Hence, one 
has 

ni = --(^], (D12) 

where It = Iq — la is the transmitted light intensity, and Id 
is the intensity of the background light recorded on the CCD 
camera without the imaging beam. Polarization effects and 
the atomic manifold level- structure of the optical transition 
used in the imaging process can bring the absorption cross- 
section to be smaller than its theoretical value given above. 
This would lead to underestimate the number of atoms. Thus, 
we performe an accurate calibration of the absorption imaging 
efficiency f45l. For the pixel i centered at {xi^yi) as denoted 
above, the atomic density is given by 

nf^=a{^,f^-^,^-\ (D13) 

where the prefactor a and its uncertainty cTq, are given by 

a = 0.112, o-c, = 0.009, (D14) 

and fiQ^ is an offset that may vary from image to image 

in) 

(hence the index n). We estimate fi^ ' as follows. For pixels 
i far away from the centre of the image, we do not expect any 
atoms. We consider quadratic frames centred on the image 
of thickness one pixel and increasing size. We then calculate 
the average of ii^^^ for each frame and take /ig^^ as the min- 
imum over all such frames and (F denotes the set of pixels 
corresponding to the frame) 

To summarize, the best estimation of the entanglement 
bound over a set of M images is given by 



odtA 




FIG. 5: Experimental sequence to obtain samples with different tem- 
peratures in the same final potential before ramping the lattices: 
Ramps of the optical dipole trap (ODT) and the optical lattice (OL). 
The value of Po corresponds to the power of the ODT at the end of 
the evaporation. It is tuned in order to realize samples of different 
temperatures. 

and our best estimation of the systematic uncertainty is 

M 

»=i (D16) 

which evaluates to 

2 2 / \ ^ M 

\ieA ' ' * / 

\ 2 — 2 

' ' ieA 

The global uncertainty can be found adding systematic and 
statistical errors in quadrature, 

which corresponds to the error bars indicated in the main text. 
Analyzing the different contributions, we have found that the 

in) 

main error sources are related to the estimation of a and /ig ^ , 
hence the error bars actually do not decrease when considering 
a larger number of images, i.e. statistical errors are negligible. 
In other words, our entanglement estimation does not require 
a large set of absorption images in order to get small errors on 
the lower bound. 



Appendix E: Production of samples with different temperatures 

In Fig. [5] we show the experimental sequence adopted to 
obtain samples with different temperature T (the temperature 
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before ramping up the lattice) in the same harmonic potential. 
As the temperature is varied through a final evaporation per- 
formed lowering the power Pqdt of the optical dipole trap to a 
final value Pq, in general samples with different temperatures 
are obtained in different harmonic potentials. For this rea- 
son, before ramping up the lattices we adiabatically increase 
in 500 ms the ODT to a fixed power Pi . In this way all the 



samples with different temperatures are prepared in the same 
harmonic potential with cylindrical symmetry characterized 
by a radial frequency of 50 Hz and an axial frequency of 8 Hz. 
The atom number in the samples with different temperatures 
is kept constant by varying the loading time of the magnetic 
trap in order to realize samples with sinular atom number. 



